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ABSTRACT 

The low-mass condensations (LMCs) have been observed within the molec- 
ular cloud cores. In this research, we investigate the effect of isobaric thermal 
instability (TI) applied for forming these LMCs. For this purpose, at first we 
investigate the occurrence of TI in the molecular clouds. Then, for studying the 
significance of linear isobaric TI, we use a contracting axisymmetric cylindrical 
core with axial magnetic field. Consideration to cooling and heating mechanisms 
in the molecular clouds shows that including the heating due to ambipolar dif- 
fusion can lead to the occurrence of TI in a time-scale smaller than dynamical 
time-scale. Application of linear perturbation analysis shows that isobaric TI 
can take place in outer region of the molecular cloud cores. Furthermore, the 
results show that perturbations with wavelengths greater than few astronomi- 
cal units are protected from destabilization property of thermal conduction, so 
they can grow to form LMCs. Thus, the results show that the mechanism of TI 
can be used to explain the formation of LMCs as the progenitors of collapsing 
proto-stellar entities, brown dwarfs, or proto-planets. 

Subject headings: ISM: clouds - ISM: evolution - Hydrodynamics - ISM: mag- 
netic fields - diffusion - stars: formation. 



1. Introduction 



The molecular gases have a hierarchical structure which extends from the scale of the 
cloud down to much smaller masses of unbound substructures. The terminology for structure 



- 2 - 



of molecular clouds is not fixed. The over-dense regions within giant molecular clouds are 
termed clumps. The massive clumps are parenting clouds for the clusters of proto-stellar 
cores. These dense cores are expected to undergo gravitational collapse evolving towards 
singular points and designated as proto-stars (e.g., Stahler and Palla 2004). Two important 
informations intended for the molecular cloud cores are their shapes and density-profiles. 
Unfortunately, the core shape cannot be easily concluded from the apparent observations of 
the-plane-of-sky because it is impossible to be directly de-projected. Although, some works 
indicate a preference in prolate cores (e.g., Ryden 1996), but there are also some studies that 
favour to a strong preference in oblate cores of finite thickness (e.g., Tassis 2007). But about 
the density-profile in the molecular cloud cores, the observations show the systematically 
certain prominent common feature: a close-to-constant density over the central region (r < 
Tin), followed by nearly a power-law decline over large distances (r > rj„) (e.g., Hung, Lai 
and Yan 2010). According to these observational results, here, wc turn our attention to the 
molecular cloud cores with cylindrical shape and appropriate decreasing density-profile. 

The observational surveys of the interior of the molecular clouds have been increased 
enormously in the last decade, thanks to the increase of resolution provided by new millimeter 
and sub-millimeter interferometers, and also they are due to the systematic combination of 
observations of dust and molecular tracers (e.g., Bergin and Tafalla 2007). In this way, 
many embedded condensed objects within each star- forming core have been revealed (e.g., 
Pirogov and Zinchenko 2008) which called low-mass condensations (LMCs). They may be 
converted to the star-forming gravitationally unstable proto-stellar cores, brown dwarfs, or 
proto-planets. Historically, we can refer to the work of Langer et al. (1995) who observed 
LMCs in the core D of Taiuus Molecular Cloud 1 in the regime of 0.007 — 0.021pc and 
0.01 — O.15M0. For the recent observations, we can refer to the discovery of very low 
luminosity objects by the Spitzer 'From Molecular Cores to Planet-Forming Discs' (c2d) 
project (e.g. Lee et al. 2009), or the results gained by Launhardt et al. (2010) in which 
they found that at least two thirds of 32 isolated star-forming cores, which were studied 
there, are the evidence of forming multiple stars. Clearly, these observations can also be 
used as a witness for existence of LMCs in the parent core before conversion to the star- 
forming entities. Although, due to limitations in resolution and sensitivity of the available 
observations (Tafalla 2008), our understanding about substructures within the cores is still 
incomplete, but this is well established that the LMCs are ubiquitous within the molecular 
cloud cores. 

The LMCs move in the parent core, thus, they may be suitable places for assembly 
the dusts and other molecules; they also may be merged with each other (e.g., Nejad- 
Asghar 2010). This topic is beyond the scope of this paper. Here, we want to investigate the 
theoretical aspects explaining the formation of LMCs within the cores. Since the molecular 
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cloud cores are characterized by subsonic levels of internal turbulence (e.g. Myers 1983, 
Goodman et al. 1998, Caselli ct al. 2002), small rotational velocity gradients (e.g., Goodman 
et al. 1993, Caselli et al. 2002), and subsonic infall motions (e.g., Tafalla ct al. 1998, Lee, 
Myers and Tafalla 2001), the turbulence alone cannot be responsible for the formation of 
LMCs. Furthermore, the associated thermal pressure in the cores dominates the turbulent 
component by a factor of several (e.g., Tafalla et al. 2004). As another mechanism for forma- 
tion of LMCs, we can refer to the work of van Loo, Falle and Hartquist (2007) who examined 
the effect of MHD waves on dense cores. They found that short-wavelength waves can play 
an important role in the generation of LMCs within a core without breaking it. Here, we can 
encourage ourselves to authenticate the suggestion that the thermal instabihty (TI) may be 
a possible favorable scenario for formation of LMCs. This idea can also be supported by the 
probes in temperature of the molecular cores in which the observations show there may be 
a deviation from uniform isothermal case (e.g., Harju et al. 2008, Friesen et al. 2009). 

After the pioneer work of Field (1965), there are many works on development of TI 
in the interstellar medium, in both numerical and analytical approaches. The main idea 
is that the rapid growth of TI leads to a strong density imbalance between the cold dense 
inhomogeneity and its low-density environment. In the molecular clouds, Gilden (1984) 
calculated the cooling function to show that molecular gas may be thermally unstable in 
environment where CO cooling dominates. The numerical calculations of Falle, Ager and 
Hartquist (2006) give the idea that TI may have an important role in formation of internal 
substructures in the molecular clouds. Fukue and Kamaya (2007) considered the effect of 
ion-neutral friction of two fluids on the growth of TI. They found the TI of the weakly ionized 
plasma in the magnetic field, even at the small length scales, could be grown up. Nejad- 
Asghar (2007) applied this idea for a one-dimensional self-gravitating magnetized molecular 
slab including the frictional heating due to ion-neutral drift. He found that the heating of 
ambipolar diffusion in outer regions of the slab is more significant than the average heating 
rates of cosmic rays and turbulent motions, thus, isobaric TI could take place in this area. 
Nejad-Asghar and Molteni (2008) used the two-fluid smoothed particle hydrodynamics to 
simulate a partially ionized one-dimensional cloud. Their gained results indicated that the 
TI can insist on the occurrence of density contrast at outer parts of the molecular cloud 
cores. Also, the two-dimensional simulations of Nejad-Asghar and Soltani (2009) confirmed 
this idea that heating of ambipolar diffusion in the molecular cloud cores may lead to TI 
and formation of LMCs. 

In this paper, we try to make an extension to this idea that TI may be responsible 
for the formation of LMCs within the molecular cloud cores. For this purpose, we consider 
a cylindrical contracting cloud as a background unperturbed arrange, with a power-law 
drop in its density-profile. In section 2, the net cooling function in the molecular clouds is 
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calculated, and happening of TI is investigated too. In section 3, the linear perturbation 
analysis with Fourier decomposition of space, is applied on the contracting cylindrical core 
to attain the criterion of isobaric TI. The effect of thermal conduction on TI stabilization is 
also investigated in Section 3. Finally, Section 4 is devoted to a summary and conclusions. 



Two key parameters to examine the happening of TI are the net cooling function and 
its time-scale. Determination of the cooling rate for an optically thick, dusty molecular 
medium is a complex non-LTE radiative transfer problem. Reviewed article by Dalgarno 
and McCray (1972) presents some of the earlier estimates of cooling rate by excitation the 
rotational levels of diatomic molecules. Since at that time, numerous authors had discussed 
various aspects of this problem. For example. Goldsmith and Langer (1978) analyzed in 
detail the cooling produced by line emission from a variety of molecular and atomic species 
at temperatures of 10, 20, and 40K and for II2 densities in the range of 10® < n{B.2) < 
lO^^m"^. Hollenbach and McKee (1979) evaluated the cooling of most coolant molecules by 
using the escape probability approximation which is equivalent to assumption of constant 
source function. The more comprehensive study of radiative cooling rate and emission-line 
luminosity in dense molecular clouds carried out by Neufeld, Lepp and Melnick (1995). 
They considered the radiate cooling of fully shielded molecular astrophysical gas over a wide 
range of temperatures (10 < T < 2500 K) and H2 densities (10^ < n{B.2) < lO^^m"^). 
Goldsmith (2001) computed the cooling effects of molecular depletion from the gas phase on 
grain surfaces in dark clouds. He parameterized the cooling function as 



and gave the value of parameters to different depletion runs. This equation was used by 
Nejad-Asghar (2007) to investigate the formation of fluctuations in a one-dimensional molec- 
ular cloud layer. 

In this paper, we use the cooling function based on the work of Neufeld, Lepp and 
Melnick (1995) which allows us to include cooling from potentially important coolants of 
five molecules and two atomic species: CO, H2, H2O, O2, HCl, C, and O. The results 
presented in figures 3a — 3c? of Neufeld, Lepp and Melnick (1995) are convenient to do rough 
parametrization like the equation ([T]), specially for temperatures between lOK and 250K. 
The effects of chemistry upon the cooling rate are dramatically illustrated by the fraction of 
the cooling rate attributed to H2O, which rises rapidly above temperature ~ 300 K due to 
an increase in water abundance. Thus, we focus our attention to the molecular clouds with 
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The occurrence of TI 




(1) 
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temperature less than 200 K. For this temperature range, the values of parameters A(„) and 
/3(„) are given in the Fig. [H 

As we know, the thermal process may be important if the dynamical time-scale is in 
excess of the cooling time-scale. We consider the contraction time-scale of a cylindrical 
molecular cloud as a multiple rj > 1 of the free-fall time-scale. 



, Stt 



32Gmn 

3.5 X 10^ year, (2) 



n 

where the latter is written for the spherical uniform density distribution. The cooling time- 
scale tc = 3kT/2mA and the free-fall time-scale ([2]) are shown in Fig. [2j In fast contraction 
(t] ~ 1), TI is important in small densities (e.g., less than lO^^m"^ for T = lOOK), while 
in slow contraction (77 >> 1) in which the cooling time-scale is much smaller than the 
contraction time-scale, the importance of TI as a trigger mechanism is much evident. 

There are several different heating mechanisms in models of interstellar matters. Since 
the ultraviolet photons are mostly screened out in dense molecular clouds, heating by col- 
lisional de-excitation of H2 molecules after radiative excitation of Lyman bands, photoe- 
mission from grains, radiative dissociation of if 2, and by chemical reactions is not im- 
portant. In addition, because of the small neutral hydrogen abundance, heating due to 
ejection of newly formed H2 molecules from grain surfaces is negligible. The heating due 
to cosmic rays with sufficient energies (~ lOOMeV) to penetrate dense clouds is com- 
monly about TcR = 2.5 x 10~^J.kg~^.s~^, with assumption of an ionization rate per H2 
molecule of 2 x 10~^^s~^ and a mean energy gains per ionization of 19eV (e.g., Glassgold 
and Langer 1973). Following Black (1987) the turbulence dissipation heating rate can be 
estimated as 

i-TR ~ 

mn 

- 1.6 X 10-«(-^)3(i^) J.kg-^s-\ (3) 
1 km.s t 

where Vturb is the turbulent velocity and I is the eddy scale. With vturb ~ 1 km.s~^ and 
/ ~ 1 pc, we obtain Ttr ~ 1-6 x 10~^J.kg~^.s~^, comparable to half of the heating rate 
of cosmic rays. In this way, we collect the values of these two heating rates to obtain 
~ 4.1 X 10^® J.kg^^.s^^. Another important heating mechanism of self-gravitating contracting 
cloud is the heating produced by gravitational compression work. An estimation for this 
heating can be derived directly from the rate of compression work per particle, pd{n'^)/dt, 
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and is given by 

Tgr ~ ^ 

where we have taken dn/dt ^ n/tcn appropriate for a uniform contracting spherical cloud 
with contraction time-scale ten = V^ff where >> 1 is for slow contraction, and ~ 1 is for 
fast one (i.e., free-fall case). 

The dissipation of magnetic energy would be considered as another heating mechanism, 
if this energy is not simply radiated away by atoms, molecules, and grains. The major field 
dissipation mechanism in the dense clouds is almost certainly ambipolar diffusion, which was 
examined by Scalo (1977) for density dependency of magnetic field in a fragmenting molecular 
cloud. Padoan, Zweibel and Nordlund (2000) presented calculations of frictional heating by 
ion-neutral drift in three-dimensional simulations of turbulent, magnetized molecular clouds. 
They show that average value of ambipolar drift heating rate can be significantly larger 
than the average of cosmic-ray heating rate. In addition, Nejad-Asghar (2007) considered 
a molecular slab under the assumption of quasi-magnetohydrostatic equilibrium, concluded 
that ambipolar drift heating is inversely proportional to density and its value in some regions 
of the slab can be significantly larger than the average heating rates of cosmic rays and the 
dissipating turbulent motions. To gain an insight on the heating rate of ion-neutral friction, 
we assume that the pressure and gravitational force on the charged fluid component are 
insignificant compared to the Lorentz force because of low ionization fraction, thus, the drift 
velocity Vd is inversely proportional to the density and directly proportional to the gradient 
of magnetic pressure (see equation [8] below). Here we choose, in a general form, Vd oc Kp~'' 
where k = A(i?^/2/io)/Ax is the change of magnetic pressure in length-scale Ax, and the 
value of b may be approximated near 3/2 (this value is from the assumption of ionization 
equilibrium with the ion density being power law of the neutral density with power 1/2). 
Nejad-Asghar (2007) adopted the value of b in the range between 0.5 and 2.0 to examine the 
isobaric TI in the regions of a self-gravitating molecular slab. The heating due to ambipolar 
diffusion is 

^AD = = 7ADep ' (5) 

p 

where = jAD^P^^'^^d is the drag force per unit volume exerted on the neutrals by ions, 
lAD ~ 3.5 X 10^°m^.kg^^.s^^ is the collision drag in molecular clouds, and we used the 

1/2 

relation pi = epn between ion and neutral densities in ionization equilibrium state with 
e ~ 9.5 X 10^^^kg^/^.m~^/^ (Shu 1992). Substituting the appropriate values for typical 
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molecular cloud cores, the ambipolar diffusion heating rate is as follows 

/ K \ 2 / -n \ -26+0.5 

where the general form Vd = i^p~^ Iiad^ is used according to equation ([8]) below. The total 
heating and cooling rates of a magnetized slow contracting molecular cloud core with typical 
value = 10 are shown in Fig. |3]for different values of parameters 6, k and temperature T. 

Now, we can turn our attention to the occurrence of TI in the molecular clouds. The 
standard criterion applicable to the interstellar gas, in the isobaric instability condition, can 
be written as 

- {^)% < 0, (7) 

where = A — F is the net cooling function, Qp = {dQ/dp)T and Qt = {d^/dT)p are 
evaluated in equilibrium state (see, e.g., the review by Balbus 1995). When the density 
along the net cooling curve increases so that a density fluctuation is formed, the above 
criterion states it will tend to be amplified because the mentioned fluctuation has a larger 
cooling than its heating rate, thus, the surrounding gas will compress it further. Fig. [3] 
suggests that without heating of ambipolar diffusion, there would not be any TI while with 
it, the instability criterion ([7]) can worthily be satisfied specially in attenuated region of a 
molecular cloud. Furthermore, this scenario can also be shown by the pressure-density plane 
which is an elegant method to describe the isobaric instability modes in thermal equilibrium. 
In this method, an equilibrium state is specified by the intersection of a thermal equilibrium 
curve and a constant pressure line. In this way, any perturbation from equilibrium state leads 
to split the gas in two phases, as the individual gas element transfer into related phase in 
correspondence with its initial sign of fluctuation (i.e., density increase or density decrease). 
Pressure-density diagram of thermal equilibrium and instability criterion are depicted in 
Fig. m This figure shows that including the heating rate due to ambipolar diffusion can 
produce a two-phase medium in outer regions (low density) of a typical molecular cloud 
core. 



3. Linear perturbation analysis 

For investigation the importance of TI in outer region of a molecular cloud core, a 
long axisymmetric cylindrical geometry is assumed to extend in z-direction with an axial 
magnetic field B = -Bq^, which directly coupled only to the charged particles. 



- 8 - 



3.1. Contracting cylindrical core 

Principally, the ion velocity Vj and the neutral velocity v„ in the molecular clouds, should 
be determined by solving separate fluid equations of these species, include their coupling by 
collision processes. However, in the time-scale considered here (see. Fig. [2]), two fluids of ion 
and neutral are approximately coupled together with a drift velocity given by 

Vd = v, -v„^ L^(VxB)xB, (8) 

which is obtained by assuming that the pressure and gravitational forces on the charged 
fluid component are negligible compared to the Lorentz force because of the low ionization 
fraction. Here, in a good approximation we choose p = Pn + Pi ^ Pn- In this way, we can 
use the basic equations as were given by Shu (1992): 

^ + V ■ Vp + pV ■ V = 0, (9) 

5v B 

pj- + p(v ■ V)v + V(p + — ) - (B ■ V)— - pg = 0, (10) 
ot zpo po 

^ ^ + ^ V ■ + ^pV ■ V + - V ■ (KVT) = , (11) 
V ■ g = -iTrGp, (12) 

p--pT = 0, (13) 
P 

where K ^ 2.16 x 10-^T^/^J.s-\K'^ .m ^ is the thermal conduction coefficient in molecular 
clouds (Lang 1986), and other variables and parameters have their usual meanings. The 
energy equation (fTTj) contains a source term — |pV ■ v, which describes the work is done by 
expansion or contraction of the medium. Thus, we note that the self-gravitating heating 
term (jlj) must be excluded from the net cooling function Q. The time evolution of magnetic 
field itself may then be obtained from the approximation that it freezes only in the plasma 
of ions and electrons. In this way, we obtain a nonlinear diffusion equation as follows 

— + (v ■ V)B + B(V ■ v) - (B ■ V)v = V X (v^ x B). (14) 

If the right-hand side of the magnetic induction equation (fT4l) is negligible, the field is well 
coupled to the whole fiuid. 

As a basis for small-perturbation analysis, we consider an axisymmetric cylindrical 
background that is contracting in a quasi-hydrostatic balance between the self-gravity and the 
repulsive forces (i.e., thermal and magnetic pressure forces). The background quantities will 
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be denoted with the subscript "0" . We consider a similarity solution for the bulk unperturbed 
fluid so that its velocity field depends linearly on the axial distance r as V(r,t) = ^^/^r, 
where is a non-dimensional parameter which presents axial contraction in the co-moving 
Lagrangian coordinate (s < 0). The basic equations fE])-f lT^ lead to 



PO(r,t) = Pc 



(15) 



Tr 



0(r,i) 



0(r-,i) 



Bcfi- 



(16) 



(17) 



where pc, Tc, and Be are respectively initial central density, temperature, and magnetic field 
at the inner radius r^n, /(;r-) is only a function of radius, and follows the equation 



1 - k-) 

where to is a time-scale free parameter (t < to)- The equations (|T5l) -( fT7j) are accurate for 
r > Tin, and we assume the density is homogenous for r < rj„. For details, the reader is 
referred to the appendix. At the end of this contraction process in which t is equal to to? 
the minimum value of Sf^t) takes place, thus the maximum effect of this contraction is an 
approximately tenfold over density in its profile. 



3.2. Linear perturbation 

For obtaining a linearized system of equations, we split each variable into unperturbed 
and perturbed components, the latter is indicated by subscript "1". Then, we carry out a 
spatial Fourier analysis with components proportional to exp(ifcr) where k is the component 
of perturbation wave-vector perpendicular to the magnetic field. Simplifying the linearized 
equations fl5]l- ffT^ by repeated use of the unperturbed background equations ffT^ - ffTS]l . we 
obtain 



^ + [tkr + 1) -pi +(tkr-^]^v,= 0, (19) 
at s \ 6 ) r 

dvi (s A'KGpor\ s fikBo , 1 dBo\ , . 

Po~J7 + -'^ - ^0 + -J — — Pi + -Povi + ikpi + + — LBi = 0, (20) 

at \s ikr + I J s \ Po Po or J 
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dBi 



where 



S 1 diWr}) 

{ikr + 1) - + ikvd + — hr^ 
s r or 



Pi 



Bi + ikBo + 



dBo 
dr 



1 ^ / O N 

Vi + -TT- [rBoVdi) 
r or 



—poTi + —TqPi, 



Vdl 



3/2 
PolADePo 



ikBn 



dB, 
dr 



SBodBo 
2 po dr 



-Pi 



0, (22) 
(23) 

(24) 



The perturbed drift velocity due to the perturbed magnetic field changes the ambipolar 
diffusion heating rate. In fact, the net cooling function Q in equation ( 12T]) must contain also 
the ambipolar heating due to perturbed magnetic field. However, because we have expressed 
the ambipolar heating in equation ([6]) by parameters, in numerical calculations for drawing 
figures, we can involve this effect by the amount of parameters. Thus, for simplicity, we 
neglect the explicit calculations of the changes of ambipolar heating due to perturbed drift 
velocity. 

Since the equilibrium is time dependent, the normal modes of the system are time 
dependent too, thus, we must apply some approximation techniques namely WKB approxi- 
mation (e.g.. Bora and Baruah 2008) to gain valuable insight into the nature of problem; this 
analysis may be a challenging task in a subsequent research. Here, we consider the isobaric 
TI, which is more realistic phenomena in the interstellar gases. Gathering the equations fl23|) 
and fl2T|) with equation f|T9|) . in isobaric case {pi = 0), leads to an exponential growth as 
follows 



pi = pi{t=o) exp 



00 {r, t')dt' 



where u{r,t') is a complex function with real and imaginary components given by 



• 1 I 35p 2 



, k' 

kx 



Co, 



and 



^[uir,t')] 



s kr 

S ^1+(|)2A;V2 



kr - ko + 



k 



K 



Co, 



respectively, where the notation of the pioneer work of Field (1965) is used as follows 



^ _2 pQ^T 
^ 3 -R Co 



^ 2 p poQp 
^"3i?ToCo' 



k 



K 



3 R PqCq 
2p Ko'' 



(25) 



(26) 



(27) 



(28) 



where Cq is the speed of sound in the unperturbed medium. The Field's isobaric instability 
criterion [kx < kp — k'^/kx) can be clearly revived via the equation f l26|) by considering a 
stationary cloud (s = 0). 
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Manipulating the equation fl26|) leads to the isobaric TI criterion in a contracting ax- 
isymmetric cylindrical cloud core as follows 

lT{r, t) - Ip{r, t) + er^lKir, t) + i±i|J!^J,(t) < 0, (29) 



l + ||A;2r2 



where 



2 u r' " 
lT(r,t) = Sir d{t'M, 



3 -R Jo ^0 

iK[r,t) = -- —— d{t/to) = -- — r — r(i-or)' 



SRJq r^po ?> R pcTin \rin J to 3 to 

m ^ / - d{t'/to) = -ln(l---), (30) 



Jo c 2 3 to' 

where the Ik and /c are explicitly evaluated, using the background equations (n3]) - (IT5]) . and 
= 2.16 X lO-^Ti/^J.s-^K-i .m ^ is the thermal conduction coefficient in the axis of the 
cylinder. 

There are three notable remarks for Ik and Ic in the equation ( !29|) : (I) Since Ik has a 
positive value and Ic has a negative value, their effects on instability criterion are reciprocal. 
Several studies have demonstrated that thermal conduction can stabilize and even erase 
a density perturbation (e.g., Burkert and Lin 2000). This implies that in the opposite of 
thermal conduction, the contraction can fortify TI because the energy is transported inward 
by the contraction process. (II) During the growth of density perturbations in isobaric 
regime, the resulting temperature gradient induces the conductive heating in fluctuations. 
The temperature gradient is greater at the distant region of the axis because in isobaric 
perturbation (Ti = — — pi), po is a decreasing function of the axial radius r. Thus, the 
larger temperature gradient at outer regions leads to amplify the Ik for more suppression of 
perturbations. (Ill) Consideration of time dependencies of Ik and Ic shows that increasing 
of \Ic\ is faster than Ik for t > 0.5to- Thus, in general the formation of LMCs via isobaric 
TI may be suppressed by the thermal conduction process as long as it overcomes on the 
contraction effect. 

Scrutiny in the isobaric TI criterion (12^ requires to apply some numerical values from 
the typical molecular cloud cores. Here, we consider a typical core with a dimension Tout = 
O.lpc and central density Uc = 4.6 x lO^^m"'^, which its total mass is approximately close to 
the one solar mass. The central uniform region is assumed to spread out to rj„ = 0.002AU. 
Thus, the density in the envelope will be nout = lO^m"'^, which is obtained from initial 



density-profile f|T5l) . Choosing Tc = 15K and b = 1.5, provides a suitable situation for 
happening the TI as depicted in Fig. HJ The thermally unstable region of a contracting 
axisymmetric cylindrical core, according to criterion f l2^ . is shown in Fig. [S] for t = to 
t = to, where it was assumed that to = IMyr. In this figure, the effect of thermal conduction 
that is true in the limit of very long wavelength (fc^r^ << 1) is ignored. We can see from the 
shade region of Fig. [5] when the thermally unstable region of the cloud shifts to the outer 
parts. If we consider an inner marginal part of the core which is thermally unstable (for 
example r = 10^'^rj„), it may be converted to a LMC via TI process. But, by the time, this 
part will be stable as is shown in Fig. [5], thus, this preformed LMC will be self-perpetuated. 

As we know, the medium may be stabilized by allowing the thermal conduction, and 
its significance depends on the perturbation wavelength which is known as Field's length 
(Field 1965). For this purpose, the area of unstable region in the r — t plane of Fig. |5]is de- 
picted in Fig.[6]for various values of /c^r^ (i.e., different perturbation wavelengths). According 
to this figure, we can realize two critical wavelengths that are kdV ~ 10^'^ and kc2r ~ 10^. If 
we consider the maximum value of radius as r ~ 10^'^rj„, the first critical wavelength will be 
Aci ~ lO^Tin- Perturbations with wavelengths greater than the mentioned critical value, in 
outer region of the molecular cores are thermally unstable and may be considered as suitable 
places for the formation of LMCs. On the other hand, using the minimum value of the radius 
of unstable region, r ~ lO^rj^, the second critical wavelength will be Ac2 ~ lO^rj^. Pertur- 
bations with wavelengths between Ad and Ac2 are partially affected by thermal conduction 
while perturbations with wavelengths less than Ac2 will entirely be suppressed. For kc2 we 
used the maximum allowed r in the unstable region and for kd its minimum was used. These 
choices are because of obtaining the minimum value of wavelength for Ac2 and the maximum 
value of it for Ad. 

It is interesting to obtain the growth time-scale of the TI as follows 

^ j:nuir,t')]dt' CO 

Ingrowth ~ T — TJT' [-'T ~ p ~r 

t t/to 

where is written for the limiting case fc^r^ << 1. The result is approximately about 0.2Myr 
that is much smaller than Jeans and magnetic instability time-scales (e.g., Fiege and Pu- 
dritz 2000, Basu, Ciolek and Wurster 2009), as we also expected from the Fig. |2l Thus, the 
TI can justify the formation of LMCs from all perturbations greater than Ac2 in outer region 
(low density) of the molecular cloud cores. 



(31) 
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4. Summary and conclusions 

In this paper, we attempted to represent importance of TI as a trigger process for con- 
figuration the hierarchical structure of molecular clouds, specially formation of the observed 
LMCs in the cloud cores. The cores are usually well thought-out as cool molecular gases 
in which the temperature and density are not completely uniform. Here, we considered the 
cooling rate of the molecular cloud (Fig. [1]) and various mechanisms of heating rates. In- 
vestigation of the net cooling rate in the molecular clouds shows that the consideration of 
the heating due to ion-neutral drift can lead to the thermally unstable gas (Figs. [3]and|l]). 
As we know, the significance of a physical process might be regarded by its time-scale. For 
this purpose, we compared the cooling time-scale and the contracting time-scale of the core, 
and we concluded that in fast contraction (free-fall), the TI is important in small densities, 
while in slow contraction, the cooling time-scale is even much shorter than the contraction 
time-scale (Fig. |2]). Thus, the TI can be considered as a trigger mechanism to imprint the 
growth of inhomogeneities within the density and temperature of the molecular cloud cores. 

For investigation the importance of the TI and finding where it will take place, we 
applied the linear perturbation analysis on a contracting axisymmetric cylindrical cloud 
core. Here, we turned our attention to the isobaric TI. The obtained instability criterion 
(129]) demonstrates that in opposite to the thermal conduction, which can suppress the hnear 
density perturbations in the medium, the contraction can fortify it because the energy is 
transported inward. Furthermore, the large gradient of temperature at outer regions of 
the core leads to more suppression of perturbations therein via the thermal conduction 
process. Another remark is that the formation of condensations via isobaric TI may in 
general be suppressed by the thermal conduction process as long as it overcomes on the 
contraction effect. Applying the numerical values of a typical molecular cloud core gave us 
some insight about the places where instability may be arisen. Fig. [5] depicts the thermally 
unstable regions of a contracting axisymmetric cylindrical core in the limitation of very long 
wavelength (i.e., the effect of thermal conduction is ignorable). According to this figure, by 
the time, unstable region of the cloud shifts to the outer region, the preformed LMCs can 
be self-perpetuated. 

The medium would be stabilized by consideration of thermal conduction which its sig- 
nificance depends on perturbation wavelength (Fig. [H]). The thermal conduction can entirely 
suppress the instability of perturbations with wavelengths less than Ac2 ~ 0.2AU, while 
its effect is completely ignorable at wavelengths greater than Ad ~ 200AU. Thus, LMCs 
more massive than mH2^^li ~ lO^^kg can grow via TI without being undergone the ther- 
mal conduction, and those smaller than mH^n\l2 ~ lO^^kg are completely disrupted by the 
thermal conduction process (n is chosen approximately lO^'^m"^ where TI occurs). Pertur- 
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bations with having wavelengths between Ad and Ac2 are partially affected by the thermal 
conduction. The relative mass which is contained within the thermally unstable region of the 

cylindrical core is approximately ^^""""2 = 1 — ( ^-^^^ ) ^ 0.99, where rmm ~ 10^'^^rj„ 

T^I'max \ '''max J 

and Tfnax ~ 10^'^^rj„ are chosen from Fig. O Thus, since approximately 99% of the core 
mass is contained within the thermally unstable region, its mass is sufficient to explain the 
formation of LMCs. Although, the TI may justify the formation of LMCs more massive 
than IO'^Mq, however, the puzzle is highly incomplete. We might investigate the behavior 
and evolution of these LMCs because they move and accumulate ubiquitous dust and other 
molecules. Furthermore, they may collide and merge to form the larger LMCs. Finally, they 
may be converted to collapsing proto-stellar entities, brown dwarfs, or proto-planets. 
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A. Background evolution 



In order to solve the set of equations fl9l)- f|T4l) . we look for a similarity contraction using 
the CO- moving system: r = S(t)X, where x is the Lagrangian coordinate and S(t) is the 
contraction parameter. In this way, the velocity field is given by V{r,t) = ^^^r. From the 
equation of continuity (|9]) , we have 

dpo , s dpo 



_l — f. 
dt s dr 



+ 2-P0 

s 



0. 



Using the separation of variables to solve the equation flAip . we obtain 

PO(r,t) = 



/ f \ 


' 1 ' 


\rinj 





(Al) 



(A2) 



where pc is the initial central density and the constraint < a < 2 fulfills the decreasing 
of density versus the cylindrical radius and its increasing with time. The equation (IA2I) is 
accurate for r > rj„, and we assume the density is homogenous (i.e. a = 2) for r < rj„. 



Substituting the equation (IA2p into the poisson equation 

-AnGpo, 



Id 
r or 



(A3) 
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we obtain, for r > r,-- 



90{r,t) 





' 1 ' 


\rinj 





(A4) 



where Qc = ATrGpcVin/a is the initial gravitational acceleration at for a long homogenous 
cylinder, and we set a = 2 for r < rin, as before. We assume that the temperature of 
the background is homogenous and it depends only on the time so that p(^r,t) = —P{r,t)T{t)- 
From the energy equation ffTTl) with the On-The-Spot equilibrium in which the energy-gain 
is locally equal to the energy- loss at that spot (i.e., ~ 0), we have 

3 1 dTo 3 1 dpQ 



2 To dt 2po dt 



3 s j^9po _^ 
2 s po dr s 



0, 



which gives 



0(r,t) 



(A5) 



(A6) 



where Tc is the initial homogenous temperature. 



Separating the time-variable from the equations of momentum (fTOl) and magnetic in- 
duction ( !T4l) by desiring the magnetic field as 



5, 



0{r,t) 



Bcfi- 



Ht) 



(A7) 



where is the initial central magnetic field strength and /(— ) is a function of radius, we 

'1^ in 

obtain the constraint a = 6 = | and the contraction parameter as 



'-to-' 



2(1 



^ 3 

2 



(A8) 



where to is a time-scale free parameter {0 < t < to). It is worth notice that the equations 
flA2l) . flA6p . and flA7p show explicitly that the ratio of magnetic pressure to gas pressure 
will be constant in duration of contraction. Clearly, this special case is emerged from our 
linear similarity assumption of the velocity field in the contracting cylinder. The separated 
r-variable of the momentum equation ( fTOl) reduces to 



2/Uo 3pcr. 



Bl 



r 



r 



+ 



^RTctl 



r 



+ 



2gct 



(A9) 



where the constant fo is chosen so that the /(^) be a real function. The equation ( ]A9[) may 
be inserted into the separated r-variable of the magnetic induction equation ( lT4l) to obtain 
via integral 



Vd = -r 

instead the approximate relation 1 



Ht) 



torf 



(-/ + r-Hr, 



(AlO) 
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Fig. 1. — The parameters of total cooling rate in the form of ^{n,T) = A(„)(T/10K)'^("' for 
molecular clouds with temperatures between 10 — 200 K. 
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Fig. 2. — The free-fall (dash) and cooling (solid) time-scales of a spherical molecular cloud 
core for three values of temperatures. 
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Fig. 3. — The total cooling (solid) and heating (dash) rates by choosing rj = 10 and k = 0.001 
for b = 1.5 and 2.0. Without heating of ambipolar diffusion (i.e., k = 0), there would not 
be any TI, while considering this heating can worthily satisfy the Field's instability criterion 
([7]) in attenuated region of a molecular cloud. 
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Fig. 4. — The Field's instability criterion (top panel) and the pressure-density diagram of 
thermal equilibrium (bottom panel) for two values of parameter b with k = 0.001, and 
without ambipolar diffusion heating (i.e., k = 0). The pressure po is evaluated at density 
lO^^m"'^ and temperature lOK. Occurrence of two-phase medium (constant pressure line) in 
the thermally unstable region is exhibited in the bottom panel as its instability criterion is 
fulfilled in the top panel. 
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Fig. 5. — The thermally unstable region (shade) of contracting axisymmetric molecular 
cloud core including the ambipolar diffusion with the parameter b = 1.5, central density 
He = 4.6 X lO^^m"^, Tin = 0.002aM, and to = IMyr. The /c^r^ is assumed to be small so that 
the effect of thermal conduction in instability criterion (p9l) is ignorable. 
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Fig. 6. — Area of the unstable region in the r — t plane of Fig. [5] versus /c^r^. The thermal 
conduction can entirely suppress the TI for a perturbation with wavelength less than Ac2, 
while its effect is completely ignorable at wavelengths greater than Ad- 



